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UNIVERSITY OF ILLINOIS 
DIGITAL COMPUTER 

LIBRARY ROUTIN E S 5 - 231 
By D. B. Gillies 

«• Natural Logarithm (D.O.I, or 3AD0I) 

Closed functional subroutine 

36 

0,1,2,3 

♦ 2"^ if y > 

8*3 to 1^«3 Millisec, depending on how much scaling is 

required . 

Enter with y > in A, link in Q. When link is obeyed, A has 

^ In y 

If y < 0, the machine will stop at 3'L on PP 02L (sexadecimal) 

with Q « y. 

Define z * 2 n y, where i < z < 1, 



and 






l+|x 
Then x — = /2 z is an algebraic identity. 

1 ' 2F X 
ln/2z.2(| + i(|) 5 + i (|)5 + ...) 

TCiis expression is replaced by an abbreviated power series 
P(x) « x + Cj x? + C 5 x 5 + Op x 7 + Cq x 9 with a maximum erroar 

of less than 10 for i < z < 1. 



Then 



j| In y - y(-nln2- § In 2 + P(x)) 



Xgp 



DAT E Fe^rua^jjfffio^T 
CODED B Y jS.jjbJfil 

checked m£)£. TKJJBi 



LOCATIOR 


ORDER 


NOTES RAGE 1 3 9 





K5 P 




J store y 


1 


if2 28L 




J plant link 




Ul F 




"~1 »*5 

J clear location for 2 ln(y) 


2 


50 IF 




y to Q for scaling 
H y 


35 F 


/ 




V 




3 


36 kL 


— j 




if y < stop 




FF 1+7F 


*-i 








4 


LO 8L 


Jl 






>• 


5 


36 10L 

00 6f 

LI 29L 


1 


\ 


is scaled y' < 2" 

— 1 6 
if so form 2 y* 

^* ■ _ 


6 


IA F 
40 F 


I 


\ 


6 -5 

and subtract ^ In 2 from 2 J In y 


7 


32 3L 


..] 




— 1 -5 

if y « 0, 2 ^ In y overflows eventually, jump to stop 


8 


22 2L 






1 loop 


01 4033F 


<r- 






1 2 y' in Q 


approximate constant 




L5 30L 








2-° 


9 


IA F 












1*0 F 






1 -*5 
^J subtract ^ In 2 from 2 ^ In y 


10 


SJ F 








H 


36 8L 






_J if y 1 < -k Junf> back. 


S5 F 






50 5L 




_J A has z ■ y 1 * scaled y| Q=*0 


12 


LO 31L 






13 


10 IF 

IA 31L 




*"-!' 




40 IF 




N(D = |(z + -±) 
\J2 


Ik 


LO 31L 






00 IF 




z - -~, small number in Q 
{2 


15 


66 IF 






S5 F 






16 


kO IF 




2(z - -±) 
1(1) - !f- * x 




TJ IF 










(* + Tr) 

>l2 





WCATIQI 


ORDER 


NOTES 


BUJE 2 




17 


1*0 2F 


| 1(2) - x 2 








50 2F 










18 


7J 32L 




C 9 x 2 








I* 33L 








19 


40 3P 




C 9 x + G 7 








50 2F 








20 


7J 3F 
L4 34L 










21 


40 3F 
50 2F 


* 


4 2 
C 9 x + c ? x + Q 






28 


7J 3F 
L4 35L 










23 


40 3F 
50 IF 




„ 6 4 p 

C 9 x + C ? x + C^ x^ + C 3 






24 


7J 3F 










25 


40 5F 
50 2F 


• 


C 9 x 7 + C ? x 5 + C^ x 3 + C^ x 






26 


L4 IF 
10 4f 




P(x) -x + C 3 x 3 + C 5 x 5 + C ? x 7 + 


C 9 x9 




27 


F4 30L 
10 IF 




^ (-In -{2 + P(x)) 






28 
29 


L4 F 

22 (link)F 


by 1 


J- ~ In 2 (- ^ In y) 
IjGbey link 

^ in 2 

-- 5±ln2 






00 F 




30 
31 


00 71,449,067,324 F 






00 F 

00 1,087, 6( 

00 F 


33,H9,889p| 








32 

* ■ 

3* 
15 


00 388,736,063,997 F 


c 9 = .ooo4 64o4 4457 

C„ « .0022 2855 8603 

c 5 * .0125 0018 5911 
C3 - .0*33 3332 94U 






00 F 

00 255,111,138 F 




00 F 

00 l,225,l63,iK)9 F 




00 F 

00 6,872,049,879 F 




00 45,312,982,352 f 
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